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Abstract. This paper investigates the thermodynamics of a large class of non- 
i-C ' Hermitian, PT-symmetric oscillators, whose energy spectrum is entirely real. The 

spectrum is estimated by second-order WKB approximation, which turns out to be 
very accurate even for small quantum numbers, and used to generate the quantum 
' partition function. Graphs showing the thermal behavior of the entropy and the specific 

^ , heat, at all regimes of temperature, are given. To obtain the corresponding classical 

partition function it turns out to be necessary in general to integrate over a complex 
"phase space" . For the wrong-sign quartic, whose equivalent Hermitian Hamiltonian is 
known exactly, it is demonstrated explicitly how this formulation arises, starting from 
the Hermitian case. 



in 
o 

1. Introduction 



PACS numbers: 03.65. Ge, 05.20.-y, 05.30.-d 



The recent surge of interest in Hamiltonians which, although non-Hermitian, 
nevertheless have a completely real spectrum began with the pioneering paper of 
Bender and Boettcher [l], which gave strong numerical and analytical evidence that 
the spectrum of the class of Hamiltonians 

H = p^- {ixf (1) 

was completely real and positive for N > 2, and attributed this reality to the (unbroken) 
PT symmetry of the Hamiltonian. Since then a large number of PT -symmetric models 
have been explored (see, e.g., [2]), and it was found that the phenomenon is rather 
general. 

The natural metric arising in such theories is not positive definite, which precludes 
a straightforward physical interpretation in terms of probability amplitudes. However, 
it turns out to be possible to construct [3] a grading operator C, which gives a positive- 
definite metric ?7+ = PC. In contrast to standard quantum mechanics, where the 
metric is the same for any Hermitian theory, here the metric is determined in each case 
by the Hamiltonian itself. In the Schrodinger wave-function formulation of quantum 
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mechanics, this metric appears exphcitly in the calculation of Green functions; however, 
in the path- integral of functional formulation its role is much more subtle [1]. 

The class of Hamiltonians ([I]) can be thought of as a continuation in the exponent 
of the potential, starting with the harmonic oscillator, = 2. A more general class of 
Hamiltonians [5l E] is obtained by continuation off the x^*^ oscillator: 

H = p^ + x^^'{zxy, (2) 

where M = 1,2,3,...; e > 0. They include the harmonic oscillator H = + 
{M = 1, e = 0), and non-Hermitian Hamiltonians such as H = + ix^ (M = 1, e = 1), 
H = p^ — x^ (M = 1, e = 2) and H = p^ — x^ (M = 2, e = 2). The eigenvalue equation 
corresponding to ^ is 

- <(X) + X^''{txYiJn{x) = Enij{x) (3) 

(throughout the text, h = = 1), where ipn{x) is required to vanish as |x| ^ oo. As 
explained in [5], Ej , when the total exponent A^ = 2M + e is greater than 4 the boundary 
condition can no longer be satisfied with x real, and one is obliged to analytically 
continue the eigenvalue equation ([3]) into the complex plane. Specifically it should 
be continued into the lower half x plane within a Stokes wedge symmetrically placed 
with respect to the imaginary axis J When this is done, the energy spectrum turns 
out to be real, discrete and positive [H El E]. A rigorous proof of this property was 
eventually constructed by Dorey et al. [7]. For reviews of the whole field of non- 
Hermitian Hamiltonians see [8]. 

In principle it is possible [9] to relate H hj a. similarity transformation to an 
equivalent Hermitian Hamiltonian h with the same energy spectrum. However, such 
a programme is difficult to implement, and an exact form for h is available only in very 
few cases. It should be stressed that, unlike ([3]), the eigenvalue equation corresponding 
to h can always be solved on the real axis. 

Possible physical applications of these unusual Hamiltonians are beginning to 
emerge. We mention here the "quantum brachistochrone" , in which the standard lower 
bound for the transition time between two states can be circumvented [lO] by a judicious 
interplay of Hermitian and non-Hermitian systems. How this could be achieved in 
practice has been recently discussed in [TT]. In analogue optical situations, where the 
refractive index plays the role of the potential, PT-symmetric systems have been shown 
to have unusual and interesting properties [12], [13] . 

Perhaps surprisingly, little attention has been paid so far to the thermodynamics 
of such systems. In the present paper we investigate the thermodynamics of the 
Hamiltonians of ([2]), which we hope will act as a template for the investigation of other 
non-Hermitian systems. As far as the quantum partition thermodynamics is concerned 
we need only calculate the energy levels, on which the partition function depends. The 
high-temperature limit then probes the classical mechanics of these systems, and we use 

I this is why the supcrficiaUy Hermitian Hamiltonians H — — x'^ and H — — are in fact 
non-Hermitian and PT-symmetric. 
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this to investigate how the classical partition function should be defined, which is not 
obvious a priori. 

The paper is organized as follows. In Sec. [2] we show that the energies En of ([3]) 
can be estimated to a high degree of accuracy, even for small values of n, using the 
WKB approximation including the sub leading contribution. In Sec. |3]we use the energy 
levels so determined to evaluate numerically the canonical partition function, which 
then yields plots illustrating the thermal behavior of the entropy and the specific heat, 
at low, intermediate and high temperatures. The plots also show how the parameters 
M and e affect these quantities. In the high-temperature limit, still using the WKB 
formula for En, we obtain a closed formula for the semiclassical partition function by 
integration over n. 

This is then used, in Sec. IH as a tool for investigating how the semiclassical partition 
function Zd associated with H in should be expressed as a "phase-space" integral. 
In general this must differ from the usual real phase-space integral, which does not 
converge. The special case H = — (the wrong-sign quartic oscillator) is ideal 
for testing these ideas, because in this case one has an explicit form for the equivalent 
isospectral Hermitian Hamiltonian h. Thus one can start by calculating the semiclassical 
partition function using the real phase space associated with h, and then transform the 
variables to find the correct integral expression for Z^i in terms of the variables appearing 
in H. Sec. [5] includes a summary and further discussion. 

2. The WKB approximation for the spectrum 

An extremely accurate approximation for the spectrum oi H in ^ can be obtained 
using the WKB method including the subleading contribution. In [H [5] the spectra 
were so calculated for M = 1. This section generalizes those results to the case M > 1. 
One begins by solving the equation x'^'^{ixY = E to determine the turning points, 

recalling that = 2M + e. It should be pointed out that x± are in the wedge mentioned 
in the previous section. For a given M, when e = the turning points are x± = ±E^^^ . 
As e increases from zero, these points migrate from the real axis towards the negative 
imaginary axis of the complex x plane. 

The leading WKB contribution is obtained by imposing the condition 

{n + l/2)7r = J^^ dx^E -x^^^ixY (5) 

over a path for which the integration is real. Choosing the path as the ray going from 
X- to 0, followed by that from to ([5]) can be recast as 

(n + l/2)7r = 2 sin(M7r /N)E^/^+^/^ dsJl - s^ , (6) 

Jo 

which yields the first factor in ([8]) below. 
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To calculate the subleading WKB contribution, instead of ([S]), one uses [HI 

{n + l/2)n = (l/2z) (f dxQ^^^ + (l/2z) (f dxQ"/ASQ^'^, (7) 

where Q[x) := x^^\ixY — E and the contour C encircles counterclockwise the rays used 
to calculate the integration in ([5]). It follows then that 

1 / dx ^"jj!;} , = —N(N - 1) sin(M7r/iV) 



;i-s^)3/2^ 



which when added to ([6]) yields (for large n) 



En 



/ r(3/2 + l/iV)v/^(n + l/2) \ 
V sin(M7r/A^)r(l + l/iV) / 



A^(A^ - 1) sin^(M7r/A^) cot(7r/iV) 
37r(n + l/2)2(Ar + 2)2 



(8) 



after identities involving products of gamma functions are used. 

The energy spectrum corresponding to the Hamiltonian p^ + lxl^ , which we consider 
for purposes of comparison, is obtained [5] by omitting the factors sin(M7r/A^) and 
sin^(M7r/A^) in ([8]) [corresponding to M = A^/2], and the results in P,[5] are reproduced 
by setting M = 1 in ([8]). Tables 1 and 2 show that ([H]) is an excellent approximation for 
both the Hermitian and non-Hermitian versions even for small quantum numbers. 



Table 1. Energy levels £"„ of the potential The successive columns are n, first- 
and second-order WKB, the numerical results from |15j . and the approximation from 





n 


WKBl 


WKB2 


Exact 


Turschner 





0.86714532 


0.98982129 


1.06036209 


1.032458 


1 


3.75191992 


3.81089637 


3.79967303 


3.785676 


10 


50.2401523 


50.2562691 


50.2562545 


50.254484 


50 


407.868707 


407.874363 


407.874363 


407.87365 


100 


1020.986417 


1020.989992 


1020.989992 


1020.989538 



Table 2. Energy levels En of the potential — x^. The successive columns are n, first- 
and second-order WKB, and the numerical results from [T]. 

n WKBl WKB2 Exact 

~0 1.37651 1.47388 1.4771 

1 5.9558 6.00261 6.0033 

2 11.769 11.8023 11.8023 

3 18.4321 18.4588 18.4590 
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3. Quantum Statistical mechanics 

In quantum statistical mechanics the partition function is given by 

oo 

Z{T) = £ e-^"/^. (9) 

71=0 

Note that this is exactly the same as it would be for a Hermitian Hamiltonian, with no 
appearance of the metric operator. As was first pointed out by Jakubsky [17], this is 
due to the cyclic property of the trace in Z = Tre~^/^. 

For the PT-symmetric oscillators we are considering, we use the second-order WKB 
approximation of ([8]) to evaluate the energies En- This is guaranteed to be accurate at 
high temperatures when the populations of states with high energies are appreciable, 
so that the relevant values of n are large. However, as we have illustrated above, 
the energies so obtained are very accurate, even for small n, so the corresponding 
approximation for Z should also be accurate at low temperatures. 

Once we have calculated Z{T) we may evaluate other thermodynamic quantities 
in the standard way from the free energy F = —TlogZ, which yields the entropy 
S = —dF/dT, the internal energy U = F + TS, and the specific heat C = dU jdT . In 
Figures [U and [2] we illustrate the thermal behavior of S and C for M = 1 and various 
values of e, and compare them with those obtained for the Hermitian potential 




Figure 1. S versus T for N =3 (top) and 10 (bottom), using tlic second-order WKB 
approximation ([5]) for the energies. In each case the solid curve (green) corresponds to 
V = |2;|^ and the dashed curve (red) to F = x^iixY'"^ ( M 1, e = - 2). 
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Figure 2. C versus T for iV=3 (top) and 10 (bottom), using the second-order WKB 
approximation ([5]) for the energies. In each case the solid curve (green) corresponds to 
V = \x\^ and the dashed curve (red) to V ^ x'^{ix)^-'^ { M = I, e = N - 2). 



Apart from the case where M = 1 and e = (the harmonic oscillator), there is 
no closed form for the summation in (jH]). However, a closed form can be obtained at 
high temperatures, where as mentioned above, ([8]) is most reliable. Working with the 
leading contribution in (IHl), since the subleading contribution is suppressed at large n, 
we obtain 

oo 

Z(r) = ^ e-'^^°(2"+^)'"''""^ (10) 

n=0 

where Eq denotes the first factor in ([8]) when n = 0, and (3 := 1/T. As [3Eq — > 0, the 
summation can be replaced by an integration, 

Z,,(T)= rrfne-(®/^)("+l/^)^"^^""^ (11) 



where G is the characteristic temperature of the oscillator, 

[sin(M7r/A^)r(l + l/iV)J ' ^ ^ 

In fact, i3Eq — > means 0/T ^ 1; i.e., one is deahng with the classical regime, as 
indicated by the change in notation from Z to Z^. A few manipulations in ( fTTI) lead to 

Z,,(T) = r(3/2 + l/iV)(^-j . (13) 
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By simple inspection, f ll3p yields the well known expression T/0 for the harmonic 
oscillator N = 2. 

Zci(T) in (fT3|) leads to the entropy 

S = [1/2 + 1/N] [log(r/e) + 1] + log r(3/2 + 1/N) (14) 

and specific heat 

C= 1/2 + 1/N, (15) 

corresponding to the (classical) thermal behavior on the right-hand part of the plots. 
The characteristic temperature B in (fT2l) gives the magnitude of the energy gap 
separating the first excited state and the ground state. As (fT2l) clearly shows, for a 
given M, G is an unbounded increasing function of A^. Since G separates the classical 
thermal behavior [ (fl^ and f|T5|) ] from the quantum thermal behavior (corresponding in 
the plots to the drop toward zero as T — > 0), one sees that the larger the higher the 
temperature up to which the quantum behavior still prevails. These features are clear 
in the graphics. 

It should be noted that 

Zci{T) = Qci{N,T) sm{M7r/N), (16) 

where 

O (NT)- + (17) 



is the semiclassical partition function corresponding to the Hamiltonian p + |x| |18j . 
As Zci and Qd are proportional to each other, they lead to the same thermodynamics 
as T ^ oo (see plots), and in particular the specific heat is the same in both cases. 
Nevertheless there is an important difference, namely, the characteristic temperature 
associated with Q^i is bounded [it is given by G in (fT2|) omitting sin(M-7r/A^)], and 
consequently their corresponding thermodynamic functions may differ significantly at 
intermediate and small temperatures. For example, if T <C G, the corresponding specific 
heat would have the usual quantum behavior expressed by an exponential decay as 
T ^ 0, 

C^f^]\-^/T, (18) 



whereas the same T might be much greater than the characteristic temperature 
associated with Qd, for which the classical behavior in (fT5l) would be observed. [Strictly 
speaking, G in ( fT8l) should be replaced by the energy gap between the first excited state 
and the ground state, but this inaccuracy does not spoil the argument.] 

4. Classical statistical mechanics 

In standard classical statistical mechanics, with a Hermitian Hamiltonian, the partition 
function is given by integrating over the real phase space, 

-PHip,x) 



]^ /■+00 r+oo 

ZcliT) = — dp dx 

ZTT J — oo J — oo 
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1 



/. 



(19) 



— oo 



for H = + V{x). 

For the non-Hermitian Hamiltonians of ([2]) it is not immediately clear what is the 
correct formulation. For sufficiently small e the integral of (fT9l) is still convergent, and 
indeed reproduces correctly the high-temperature limit of (1161) . However, for larger 
values of e, for example e > 1 in the case M = 1, the integral along the real x-axis 
diverges, and in order to obtain a convergent result the contour of integration must be 
continued into the lower half of the complex x-plane. 

The situation is similar, but not identical, to the Stokes wedges encountered in the 
quantum problem [1]. Thus, let us set x = re~*^. The condition for convergence for 
large r is that Re > 0, with V = leading to 



This leads to wedges in the complex x-plane where the integral is convergent and other, 
forbidden, wedges where it is not. The forbidden wedges rotate downwards and become 
narrower as increases, in a similar manner to Stokes wedges. 

In the case M = 1, i.e. continuing away from the harmonic oscillator, the right- 
hand forbidden wedge is centered on 6^ = (^/2)(1 — and has opening angle ir/N. 
There is a mirror image in the left-half plane. So for = 3 the lower edge of the 
forbidden wedges lies precisely on the real axis, confirming that zx^ is a limiting case. 
For N = 4, the forbidden wedges include the real axis, which is no longer viable as an 
integration contour. When = 5 the upper edge of the wedges lies on the real axis, 
and thereafter the wedges lie entirely in the lower half plane. This latter situation is 
illustrated in Fig. [3] for = 6. In this case the real axis is again a viable contour, 
but in fact it corresponds to the Hermitian potential V = \x\^, whereas the hyperbolic 
contour is a possible contour for the non-Hermitian theory: it has been pushed off the 
real axis by the continuation in e, and the two contours are separated by a forbidden 
region. There is in fact another forbidden wedge centered on the negative imaginary 
axis, which is avoided by the contour shown. 

Moreover, these integrals are not only convergent, but they correctly reproduce the 
semiclassical result of ( fT6l) . Thus, for example, in the case A^ = 6, if we integrate along 
the real axis we obtain Qci{Q,T), corresponding to the Hermitian case M = 3, e = 0, 
but if we integrate along the rays x = ± rexp(=l= iiv/S), the centres of the next allowed 
left and right wedges, we obtain precisely the extra factor cos(7r/3) = sin(7r/A^) required 
to give Zci for M = 1. 

In principle, an alternative way of determining what is the correct expression for 
the classical partition function is to start with the definition in terms of the equivalent 
Hermitian Hamiltonian and make the appropriate changes of variables. 

As mentioned in the Introduction, in quantum mechanics a non-Hermitian 
Hamiltonian H with a completely real spectrum is related to an isospectral Hermitian 
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Figure 3. The classically forbidden wedges (hatched) for = 6. The hyperbolic 
contour is a possible contour for the non-Hermitian V ^ {M = 1, e = 4), whereas 
integration along the real axis corresponds to y = \x\^ (A/ = 3, e = 0). 



Hamiltonian h by the similarity transformation [9]: 

H = p-^hp = e^^he--^^, (20) 

where the metric 77+ is written [3] in terms of the Hermitian operator Q as 77+ = e~*^. 
The operator p is therefore Hermitian and positive definite. Thus we can certainly write 
Zci as 

Z^iiT) = ^ ir dxdp e-^'^(^'f). (21) 

From (|20l) we have the condition of quasi- Hermit icity on H: 

= e-^He^ = 7]H7]-\ (22) 

Similarly any observable A, with real expectation values, must also be quasi-Hermitian, 
i.e. 

= e-'^Ae^ = r]Ar]-\ (23) 

This in turn means that A is related to a standard Hermitian counterpart a by the same 
similarity transformation: 

A = e^^ae'^^ . 

We can use this latter equation to derive two different relations between H and 
namely 

H{x,v) = e^^h{x,p)e~^^^ = h{X,P), 

(24) 

H{X\P^) = e^^h{X\P^)e~^^ = h{x,p), 
It is the second of these identities that we need here, to write 

-/3_H'(Xt, pt) 



ZcliT) = ^ jj dxdp e 
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Now make the change of variables to ^ = X\ vr = P''', to be treated here as classical 
variables. Thus, in the first instance (the Jacobian is 1), 

Zci{T) = ^ If de^^vre-^^^'-), (25) 
Ztx jjc 

where C is a contour in complex (,^, n) space determined parametrically by = X^(x,p), 
TT = P''"(x,p). That is, the initial contour is ultimately determined by the metric t]. 
Whether or not the contour can subsequently be deformed to a standard real phase 
space (^, vr) is a matter to be determined for an individual Hamiltonian. From our 
preceding discussion it seems clear that this is not possible in general because of the 
presence of intervening forbidden wedges. 

Note that the metric rj does not appear explicitly as an integration measure in 
the integral representation ( l25i) of the partition function, only through the form of the 
relationship between (^, vr) and (x, p). The issue of the role of t] in path integrals has 
been addressed and clarified in [H [19] . 

As mentioned above, there are very few cases where Q, and hence the relations 
between [x, p) and {C,, vr), are known exactly. Apart from the rather trivial case of the 
Swanson model [20] (where in fact standard phase space can be used), a nice example [21] 
within the class of Hamiltonians of ([2]) is the wrong-sign quartic V = —x^. 

A word of clarification is in order here. The original Hamiltonian in this case is 
really Hz{z,Pz) = pI — z'^, where the eigenvalue problem has to be posed on a complex 
contour in the appropriate Stokes wedges. It is only when a particular contour is chosen, 
with the parametrization z = —2i^y{l+ix) in terms of the real variable x, that we obtain 
the non-Hermitian Hamiltonian H{x,p) for which Q was found, namely 



H{x,p) = i {l + ix,p^} - - a{l + ix) 



where, in the present case, a = 16. The Q operator is 

g = -|^ + 2p, 

which results, via the first equation of (l24l) . in the equivalent Hermitian Hamiltonian 

\ , 2 

h{x,p) = 4^ ~ 2^ + "^ • 
This is still not a conventional Hamiltonian, but becomes so on taking the Fourier 
transform. 

It is easily verified that using this h{x,p) in (!2T!) correctly gives the appropriate 
result (M = 1, e = 2) in ( fT6l) . This is true whether or not one includes the linear term 
—p/2, which in fact is a quantum anomaly [22] proportional to h. 

If we now make the transformation to the variables X"^ = ^ and = tt, using the 
second equation of (Ell), we obtain 



(26) 

i^_tt(l + ze)2. 
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We now have an expression for Z^i of the form of fl25p . in which the contour is given 
parametrically by fl26|) in terms of the real variables x and p. However, it is readily 
verified that there is no obstruction to deforming the contour to the real axis, so that 
Zci can be expressed as a real phase-space integral. But H{^, vr) has precisely the same 
form as H{x,p), and we can take the reverse step to the original variables z and pz- 
That is, 

z = - 2i^y{l + ix) 

Pz=p^il + ix), (27) 

a canonical transformation from {x,p) to {z,Pz)- For any finite x, the argument of the 
variable p^ is less in modulus than 7r/4, so that it lies within the wedge including the real 
axis that guarantees convergence of the integral / dp^ exp(— Hence the Pz integral 
can be deformed to the real axis. However, for the z integration {N = 4) the contour 
cannot be so deformed because of an intervening forbidden wedge. This is in agreement 
with the results found previously when we performed the integrals over appropriately 
chosen rays, to obtain the correct semiclassical result of f|T6|) . 

5. Discussion 

To summarize, in this work we addressed the quantum and classical statistical mechanics 
of the class of non-Hermitian PT-symmetric oscillators of ([2]). These are in principle 
related by a similarity transformation to an equivalent Hermitian Hamiltonian h. 
However, for the quantum partition function Z{T) one only needs the energy levels 
of H, which we evaluated using the WKB approximation including the first subleading 
correction. The main qualitative difference from the Hermitian oscillators (M = N/2) 
with potential turns out to be that because of the factor sin(M7r/A^) in the 
denominator the characteristic temperature Q[M,N) of (fT2|) grows without limit as 
increases, so that the onset of semiclassical behavior is progressively delayed. At 
high temperatures the entropy has the same form as that for the oscillator, except 
for the different value of 0, and the specific heat is the same, depending only on A^. 

The semiclassical partition function Zci{T) was determined in the first instance as 
the high-temperature limit of Z(T). We then investigated how this result could be 
reproduced by a purely classical calculation, and found that this could be achieved only 
by extending the classical phase-space integrations into the complex plane. Specifically, 
because of the simple dependence of H on p, the p integration can remain on the 
real axis, but the x-dependence means that the x integration can only be done within 
certain allowed wedges in the complex x plane, corresponding to different values of M. 
By integrating along the rays at the centre of the allowed wedges we verified that we 
indeed reproduced the semiclassical result of (fT6!) . In order to understand this from an 
alternative perspective we considered the special case of the "upside-down" quartic, the 
one example of the class ([2]) where h is known explicitly. Starting with the conventional 
Hermitian formulation for Z^ in terms of h we implemented the similarity transformation 
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to re-express it in terms of H, showing that indeed it required a complex contour in x 
of the type we had previously found. 

It has been shown [23] that the classical Hamiltonian dynamics for systems of 
the type ([1]) and ([2]) can be formulated in a real 4-dimensional phase space, while a 
generalized canonical structure for non-Hermitian classical dynamics has recently been 
derived in [21]. It would be interesting to see how our results for the semiclassical 
partition function can be derived in those formalisms. 
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